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Abstract 

We develop a systematic semiclassical approximation scheme for quantum 



■ Hall skyrmions near filling factors v = 2 n+i > wri ich is exact in the long wave- 

00 ' 



length limit. We construct a coherent state basis for the Hilbert space of 



ON 

Chern-Simons gauge fields and composite bosons with spin. These states are 

o. 

projected to the physical gauge invariant subspace and their wavefunctions 



explicitly evaluated. The lowest Landau level (LLL) condition is shown to be 



B 

*^ • equivalent to an analyticity condition on the parameters. 

a: 

The matrix elements of physical observables between these states are 



shown to be calculable in the limit of small amplitude long wavelength density 
fluctuations. The electric charge density is shown to be proportional to the 
topological charge density if and only if the LLL condition is satisfied. 

We then show that these states themselves form a generalised coherent 
state basis, parameterised by the values of physical observables. The theory 
can therefore be written in terms of these gauge invariant bosonic fields in 
the long wavelength regime. The off diagonal matrix elements of observables 
in these coherent states are computed and shown to vanish in the long wave- 
length limit. Thus we are able to prove that the classical description of the 
skyrmion is exact in the limit of large skyrmions. 
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I. INTRODUCTION 



The combination of a low value of effective mass and a low value of the Zeeman coupling 
in GaAs makes the spin degree of freedom relevent for quantum Hall systems in this material. 
Initial calculations by Chakraborty and Zhang [[]]] showed that the low energy quasiparticles 
in such systems were spin reversed. Kane and Lee || then argued that the quasiparticles 
could be extended spin textures that are described by topological solitons of the non-linear 
sigma model (NLSM): skyrmions. Sondhi et. al. |3| estimated the energy of skyrmions using 
an effective NLSM and showed that they are energetically favoured over single spin reversed 
quasiparticles. This was corroborated by Hartree-Fock calculations by Fertig et. al. [Q. 

Experimental evidence for the quasiparticles being extended spin textures initially came 
from Knight shift measurements of Barret et. al. ||, optical magneto-absorption experiments 
of Aifer et. al. || and tilted field transport measurements of Schrieffer et. al |7|. These 
experiments found the average value of the spin of the quasiparticle to be ~ 2 — 7 depending 
on the system parameters. Subsequently, experiments have been done in systems where the 
Lande g factor is reduced by high pressures || or by engineering the material yielding 
evidence of quasiparticles with very large, (~ 20 — 30), values of spin. A negative result 



has also been reported by Kukushkin et. al. ||10|| , where optical measurements found the 
quasiparticles to have spin |, i.e. no spin textures. There has been no explanation of this 
result in the literature so far. 

In the NLSM approach, the quasiholes (particles) are assumed to be described by the 
classical solutions of the model. The electrical charge density is assumed to be equal to 
the topological charge density, thus the quasiparticles are the classical solutions in in the 
topological charge, Q top = 1, sector. Near v — 1, the NLSM energy functional and the 
relation between the topological charge density and electrical charge density has been derived 
by many workers in the LLL, long wavelength approximation P,|TT].|12| . 

What is the regime of validity of this model and is there a limit where the classical 
approximation is exact ? This is the main question we address and answer in this paper. 
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This question has previously been discussed by Girvin et. al. [|lt|, where they numerically 
compare the skyrmion energies form the NLSM, Hartree-Fock and exact diagonalisation. 
The energies match in the limit of large sized skyrmions (with size R > 15l c , l c being the 
magnetic length. This calculation suggests that the classical approximation may be exact in 
long wavelength limit. In this paper we start from the microscopic theory and analytically 
address the following questions: 

1(a). When is the topological charge proportional to the electric charge ? 
1(b). When are the corresponding densities proportional to each other ? 

2. What is the limit in which the classical approximation is exact ? 
These questions automatically throw up a third one: 

3. How can the LLL condition be imposed in the classical theory ? 



The composite boson formalism |14] is the obvious choice of the formalism to use to 



address the question of the classical limit. However, we would also like to impose the LLL 
condition, which is a condition on the quantum states. We therefore develop a coherent 
state formalism for composite bosons. 

In section |Q[ we construct a coherent state basis for the Hilbert space of the bosons and 
the Chern-Simons gauge fields. The anticommuting electron operators are then constructed 
in terms of the bosonic fields. This construction provides the explicit mapping between the 
electronic Hilbert space and operators and the gauge invariant sector of the composite boson 
theory. 

Section [HI] concentrates on the gauge invariant sector of the theory. The coherent states 
are projected on to the gauge invariant sector and their wave functions calculated. We then 
find that the wave functions thus obtained are exactly of the form written down previously 



by Ezawa ||16|| . The LLL condition is easily seen to be equivalent to an analyticity condition 
on the parameters labelling the coherent states. The coherent states are labelled by the 
values of the bosonic spinor field <fi a (x) and the Chern-Simons gauge field ai(x). The states 
projected to the gauge invariant sector depend only on gauge invariant combinations of these 
fields i.e the longitudinal part of a{x) gets related to the phase of <p(x). An inspection of 



the wavefunctions reveals another local invariance that relates the transverse part of a(x) to 
the magnitude, <f>'(x)<f>(x) of the bosonic fields. This transformation is not unitarily realised. 
Nevertheless it implies that the physical states can be labelled by a single bosonic spinor 
field which we denote by W a (x). We then show that the matrix elements of the observables 
between these projected coherent states can be computed in the limit of small amplitude, 
longwavelength density fluctuations, which we refer to as the hydrodynamic limit. We then 
discuss the relation between the charge and the topological charge densities. Our approach 



is related to previous approaches of both that of Ezawa |16[ and to that of Murthy and 
Shankar Our coherent state wavefunctions are of the same form as in reference fl6 



and our calculations are done in the same physical regime (the hydrodynamic limit) as in 



reference 15 



The issue of the classical limit is dealt with in section [TV]. We show (in the hydrodynamic 
limit), that the projected coherent states satisfy the properties of generalized coherent states 
We also show that the off diagonal matrix elements of the observables vanish in this 
limit. The theory thus admits a classical description in this limit. Further, it can be 
written in completely in terms of bosonic fields corresponding to gauge invariant physical 
observables. 



The concluding section VI summarises and discusses our results. 



II. COMPOSITE BOSON THEORY AND COHERENT STATES 

In this section, we do the usual bosonization wherein the electronic theory is written in 
terms of bosons attatched to fluxes of Chern-Simons fields. The electron operators and the 
procedure of flux attatchment is very transparent in terms of coherent states. So we first 
define the Hilbert space of the composite bosons and construct the coherent state basis for 
it. We then construct the electron operators and obtain the explicit mapping between the 
electronic Hilbert space and observables and the gauge invariant states and observables in 
the bosonic theory. 
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A. The Composite Boson Hilbert Space 



The bosonic degrees of freedom are described by the spinor field operators 
(p a {x)i fil-ix), o" = 1,2 and the Chern-Simons gauge fields, ai(x), i = 1,2. The (p op- 
erators act on the Hilbert space, Tig and satisfy the canonical commutation relations, 







The Chern-Simons gauge fields act on the Hilbert space TCcs and satisfy, 

he 



—eijd (x-y) 



K 



(1) 



(2) 



where k = 2 nhc(2n+i) • ^ we define the complex fields, a{x) and a(x) as, 

, . d2(x) + iaAx) r~K~ 

a(x) = 7= J — 

v ' v/2 \ he 

_. . a 2 (x) — iai(x) [~k~ 

aix) = -= \ — 

w v/2 She 

then the commutation relation between them is given by, 

[a(x),a(y)} = 5 2 {x-y) 



(3) 



(4) 



The full Hilbert space of the composite boson theory is the direct sum of the above two 
spaces and we denote it by, 



Hcb = 7~Lb © 7~(-cs 



(5) 



We denote the gauge invariant sector of this space by TL p h y C 'Hcb- 7~L P h y consists of the 
states which which respect the Chern-Simons Gauss law constraint, 



G(xM) 



phy 







where G(x) are the generators of gauge transformations given by, 

G(x) = /tV x a{x) — e$ (x)(p(x) 



(6) 



(7) 



We will refer to the gauge invariant observables, the operators that commute with G(x) as 
physical obsevables. 



B. Coherent State Basis 



In this section, we construct the coherent state basis for TLcb- The displacement opera- 
tors are defined to be, 

U{if) = e lj^ x ^ x ^ ~^ x ^ x ^ (8) 

where, a(x) = g^gj+Mgugj yiT ^he coherent states \a, if), parameterised by the gauge field 
a(x) and the spinor field <f(x) are then given by, 

\a,if} = U(if)D(a)\0) (9) 

where, 

a(x)\Q) = ff {x)\O) = (10) 

The states defined in equation(|9|) can be intepreted as gaussian wave packets peaked 
around the classical field configuration (a(x),if(x)). They satisfy the three standard prop- 
erties of coherent states [I7j namely, 

1. Resolution of unity: 

J T>[a,<p] \a,<p)(a,<p\ = I (11) 

where V[a,<p] = i^ M*W*)W*W*) 

2. Continuity of overlaps: 

/ I \ — f ai(x)xct2(x) — i#- f (ai(x)—a2(x)) 2 

e 5 5 x VPl( x )V2{x)-Vi{x)(p2{x)} e ~± f x \ Vl {x)-ip2{x)\ 2 

3. Values of Observables: 

(a,if\ : 0(a,a,if,if j ) : \a,<p) = 0(a,a,tp,<p ] ) (13) 
The coherent states are not gauge invariant. Under gauge transformations, 



= e t&J. a W x ™W\a-Vn,<pe-& n ) (14) 

We will now constuct a projection operator that projects any state into the gauge invari- 
ant subspace, TC P h y - Consider, 

P = — [ e ^i n{x)d(x) (15) 

where G(x) is the generator of gauge transformations, as given in equation^) and Vq = Jq, 
is the volume of the gauge group. Shifing the integration variable Q by (3 in the projection 
operator P, 

e 4-j x mG(*) P = P (16) 



Taking /3 — > 0, 



G(x)P = => G(x)P\ip)=0 (17) 



This proves that P is an operator that projects any state into H p h y - 

The above three properties (JTT] - |13D and the projection operator defined in equation ([15]) 



can be used to derive the path integral representation of the gauge invariant evolution 
operator. This is done in the Appendix [A| to obtain, 

Z = J P[a (x,t)]P[a i (x,0]^[^(^,0]e^ dtd2:rZ:(:r ' t) (18) 

where C(x,t) is the standard lagrangian of matter fields coupled to Chern-Simons gauge 
fields. This confirms the equivalence of our formalism to the standard lagrangian formalism. 



C. Bosonization 

We will now construct gauge invariant anticommuting operators that create and annihi- 
late flux carrying bosons. These operators satisfy the fermionic canonical anticommutation 
relations and can hence be used to represent the electron creation and annihilation operators 
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in Hcb- We will thus be able to map the gauge invariant sector of the composite boson 
Hilbert space, TC p h y , to the Hilbert space of the electronic system, Tt e i. The mapping is then 
used to map the observables of the electronic system to gauge invariant operators in TicB- 
We define cl(x) as 

4(z) = D(x)&(x)K(x) (19) 

We have used D(x) as short notation for D(a v x ). a v x is the classical configuration of a vortex 
with a delta function flux density at the point x. 

kV x a v x (z) = ed 2 {z-x) (20) 

D(x) therefore creates a gaussian wave packet peaked around this classical vortex configu- 
ration. When rf{x) acts on a state, ^(x) creates a bosonic particle at x and D(x) attatches 
Chern-Simons flux to it. The operator K(x) gives the Aharanov-Bohm phase corresponding 
to all the other particles already present in the state. It is defined as, 

K(x) = e**- 2 ™ 4 " 1 - 1 L e{ - X ~ z ^^^ (21) 

where 6(x) is the angle the vector, x, makes with the x-axis. 

Using the commutation relations given in equations(|l] and [|), it can be verified that the 
following canonical anti-commutation relations hold good. 

{c a (x),cl(y)} = 6 aa ,6 2 (x-y) (22) 

{c a (x),c o ,(y)} = {cl(x),cUy)} = (23) 

Hence c\(x) and c a {x) provide a representation of the electron creation and annihilation 
operators in TLcb- 

Under gauge transformations, 

<p„{x) - e^W^ix) , 0l(x) - e^^^Ux) 

Oi(x) -> Oi(x) + diQ(x) , D(x) -> e^^Dix) (24) 
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We see that c a (x) and c£.(x) are gauge invariant. 

We are now in a position to map 7i e i into 7i p h y - We map the state with number of 
electrons, \0) e i, to the vacuum state of Hob, defined in equation(|TO|), projected to TC p h y - 

\0> el -> P\0) (25) 

Since c a (x) are gauge invariant, they commute with P. Then from equation ([10]) it follows 
that, 

c a (x)P\0) = (26) 

The state with N electrons at (xi, x 2 , ....Xn) with spins (<7i, cr 2 , ctjv), |{x n , a n }jv), is then 
mapped into, 

N 

\{x n ,a n } N ) - II4„(^)P|0) (27) 

n=l 

Since the states in the RHS of equations(p5[) and ( |27| ) form a basis for 7i eZ , these equations 
specify the explicit mapping of H e i into TC P h y - 

It is now easy to map the observables as well. The density is given by, 

p(x) = cl{x)c a {x) = <pl(x)0 o (x) (28) 

The spin density is, 

S a (x) = -cKx^aa'Ca^x) = -tfl(x)T (r<T >'Pa(x) (29) 

The current density is, 

1 e 
Ji( x ) = ^( c l( x )[-^9i Ai(x)]c a (x) + h.c) 

— c 

1 c c 
= -(cpl(x)[—ihdi a,i(x) Ai(x)](p a (x) - 

— c c 

-0l(x)&(x) f aUz)G(z) + h.c) (30) 

C J z 

The last term acting on physical states is zero. Thus for matrix elements between physical 
states, we have, 



Ji( x ) = 7;0l{ x )[- i hdi — a i( x ) — A(x)]0a(x) + h.c (31) 
Similarly, the kinetic energy density, T(x) is computed to be, 

f(x) = ^-0l(x)[-ihdi - -Oiix) - -Ai(x)} 2 a (x) (32) 

III. GAUGE INVARIANT BASIS STATES 

In this section, we study the coherent states projected into Ti p h y . We show that these 
states form a basis of Ti p h y - Their wavefunctions and expectation values of observables 
are computed. The LLL condition can then be seen to be equivalent to an analyticity 
condition on the parameters. We then discuss the relation between the charge density and 
the topological charge density. Finally, we describe the parametrisation of the projected 
coherent states in terms of a single complex spinor field W a (x) discussed in the end of 
section | and derive expressions for the observables in terms of W a (x) 



A. Projected coherent states 

Consider the set of coherent states, projected to Tiphy, 

\a,if) p = P\a,ip) (33) 

Using the fact that P 2 = P and equation([TT|), we have, 

Jv[a]V[<pi]V[<p] \a,cp) pp {a,cp\ = P (34) 

P is the identity operator in Tiphy so the projected coherent states form a basis for it. 

The coherent states are not eigenstates of the number operator. Thus they have a non- 
zero overlap with states containing any number of particles. The wavefunction in the iV 
particle sector is the overlap with the states given in equation (j^), 

tp N ({xi,a,i}) = ({xi,ai} N \a,ip) p (35) 
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Using equations (|9D , Ql5|) , (|19D and (]27|), the RHS can be written as, 

i>j G ^ i=l 

e -3S5 /.(EL Si(x))x(ff(*)-vn(x)) e -^ /JsW-vnw-Etx *(*)) a ( 36 ) 



The details of this calculation and what follows is given in Appendix |B[ The Q integral in 
the RHS above is gaussian and can be done exactly. After some algebra, we obtain, 

xn^^)^^' 1,1 ^' 1,1 "'" 111 ' 11 ^^^) (37) 
1=1 

where ^({a^}) is the Laughlin wavefunction, 

MM) = IK** - ^) 2n+1 e~^ Elkl ' 2 (38) 

i>j 

a has been written as 

ai(x) = eijdjQrix) + diSl L (x) (39) 

and Qt(x) = —^fjpz. Note that when (p a (x) = constant and £lr(x) = Qt(x) =^ Vxa = B, 
the wavefunction reduces to the Laughlin wavefunction. Thus the "mean field" state is the 
Laughlin state. In this case the Q integral in equation ([3(]) is equivalent to an N vertex 
operator correlation function in a c=l conformal field theory. These wavefunctions are 
exactly of the form written down by Ezawa fl6 |. 



B. Parameterisation and LLL condition 



Apart from an overall factor that only affects the norm, the wavefunction in equation (57]) 
depends on the parameters a and cp through a spinor field W a (x) defined as, 

W (x) = ^^yU^i^-siT^-i^i)) (40) 

W a (x) and hence the wavefunction is gauge invariant (as it should be), since under gauge 
tranformations, 

11 



Qt(x) — > &t(x) 
Q L (x) -> SIl(x) + £l(x) 

^(x) - ^(x)e^ nW (41) 

a and tp have 6 real field components. The gauge invariance of the wavefunctions reduces 
the number of parameters to 5. There is another local invariance of W, i.e. 

Q T {x) -> Q T {x) + X {x) 
Q L (x) -> 

<^(x) -> ^(x)e-^ x ^ (42) 

Only the norm of the state changes under this transformation and the physical state remains 
the same. Clearly this transformation is not unitarily implemented in Tics- Nevertheless 
it reduces the number of independent real fields that parameterize the states to 4, the 
components of the spinor field W. Thus we can define the normalised projected coherent 
states, that are parameterised by W as, 

\W) = jf\a,<p) p (43) 
where, N = p (a, ip\a, ip) p , is the norm of \a, <p) p . 



From equations (p7Q and (|40| ) it is clear that the LLL condition is equivalent to the 
condition that W is analytic, 

d- z W a (x) = (44) 

Thus the LLL condition is easily implemented in this formalism as it is equivalent to an 
analyticity condition on the parameters. 

C. Observables 



We will now compute the expectation values of gauge invariant operators in the projected 
coherent states. This is given by, 
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(O) = (W\0\W) (45) 

where O is a gauge invariant observable. 

We do all our calculations the limit of W a (x) being a slowly varying function of x (over 
a length scale of l c ). As we will see, this is also the limit of small density fluctuations. We 
refer to this limit as the hydrodynamic limit. We note that this is also the limit in which 
the analytic calculations of Murthy and Shankar [[DJ are done. 

Just as in the case of the Laughlin wavefunction, the computation of M reduces to the 
computation of the partition function of a classical 2-d plasma problem. Except that here, 
the plasma density is coupled to an external field which is a function of W . In the hydro- 
dynamic limit, the partition function can be evaluated by the saddle point approximation. 
The details of the calculation are presented in Appendix 0, where we evaluate the norm to 
be, 

Af[W,Q T ] = const x e -/.^(*W»)«-*^ w -^ W} e -ft/.vnr(.)-vn r (.) 

e -8^/JM lvt ( I )»'W)+^^^)]v 2 [in(ivt W iyw)+2^n T W] ^ 

Expectation values of observables similarly reduce to the computation of expectation values 
in the plasma problem. These are also computed in the saddle point approximation in 
Appendix y. We get the density to be, 

p{x)-p = (W\p(x)-p\W) 

~ l V 2 \n(W\x)W{x)) (47) 



4vr(2n + 1) 

where p is the mean density. 

Similarly, the spin density is computed to be, 

s a {x) = (W\s a {x)\W) 
p(x 



Z\x)r a Z(x) (48) 



where we have denoted the normalised spinor by Z, 

13 



ZM - , ^ (49) 

and r a are the Pauli spin matrices. The current density and the kinetic energy density are 
also computed to yield, 

Ji(x) = (W\J t (x)\W) 

= p{x)[hLl{x) - -{Ai{x) - ai (x))] (50) 



where L'f = ^(Z^diZ — h.c) and kV x a(x) = ep( 



x 



T{x) = {W\f(x)\W) 



hu c p(x) -=- , , TTr . , (51) 



Note that the kinetic energy density is zero when W is analytic. 



D. Charge and Topological Charge Densities 

The topological charge density is given by 

q(x) = — €ijh{x) ■ dih(x) x djh{x) (52) 

where h{x) is the local direction of spin polarization, s(x) = ^p(x)n(x). In terms of Z, it is 
given by, 

q(x) = —e ij d i Z^(x)d j Z{x) (53) 



As can be seen from equations fl47|) and (|53D , the topological charge density is not necessarily 
proportional to the electrical charge density. In fact, in general, they are independent of 
each other since W'(x)W(x) and Z a [x) are independent variables. However, if the LLL 
condition is satisfied, then the analyticity of W a (x) relates the modulus and the phase of 
each component. Then W^(x)W(x) and Z a (x) are no longer independent. In fact if we use 
the analyticity condition, diW a (x) = — i€ijdjW(x), in the RHS of equation (pE7|), we get, 
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p{x) -p = - ^-— q{x) (54) 

Thus the topological charge density is proportional to the electrical charge density if and 
only if the LLL condition is satisfied. The relation fl54]) will therefore not be true in presence 
of Landau level mixing. 

When the densities are proportional, the total excess charge, Q, will of course be propor- 
tional to the total topological charge, Qtop- However, the total charges could be proportional 
without the densities being so. We will now investigate this possibility. Integrating equation 
([17]) over all space, we have, 

Q = 4 J + 1) f dx'eydj \n(W^x)W(x)) (55) 

where the contour is at infinity. If W is analytic at infinity, then the RHS of equation fl55|) 
can be written as, 

Q = — 1 I dx i \-\Z\x)d i Z{x) - d i Z\x)Z{x)) 

= -^n Qt °* (56) 

Thus if there is no Landau level mixing in the ground state, the total charge is always 
proportional to the topological charge. Note that Z and hence q(x) is well defined only if 
p(x) is non-zero everwhere. So all our considerations are true only in this case. They will 
not hold for polarized vortices where p(x) will vanish at some point. 



IV. GAUGE INVARIANT BOSONIZATION 

In the previous section, we saw that the projected coherent states are labelled by a 
spinor field W, and that the expection values of observables could be computed in the 
hydrodynamic limit in terms of W. The states can therefore be labelled by the values of the 
physical observables, the density p(x) and the normalised spinor Z a {x). In this section, we 



will show that these states themselves satisfy the generalised coherent state properties [T? 
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in Tiphy Namely, the resolution of unity and continuity of overlaps. This implies that, in 
the hydrodynamic limit, the original electronic theory can be expressed completely in terms 
of bosonic field operators corresponding to p(x) and Z a (x). Thus in this limit, the theory 
can be bosonized in a gauge invariant way with no redundant degrees of freedom. 



A. Resolution of Unity 

The fact that the identity operator in Tt p h y can be resolved in terms of the projected 
coherent states has already been shown in equation (0). Here we express this same equation 
in terms of the gauge invariant parameters. We perform the following change of variables 
in equation (|34|). 



ai(x),ip a (x) -> Q L (x),Q T (x),W a (x) (57) 

further, using equations ( pE3|) and (0), we get 

I = J T>[a,<p] \a,(p) p p (a,(p\ 

= const Jv[Q T {x)]V[Q L {x)]V[W]e-l^ nT{x W[W,Q T }\W)(W\ 

= const fv[W]g[W]\W)(W\ (58) 

where the factor e the Jacobian due to the change of variables tp — > W, 

V[W] = I| dW ^f^ (59) 

and 

G[W] = J V[tor]e- L ^ nT{x) U[W, Q T ] (60) 

The integral over Qt can be done in the hydrodynamic limit, the details are given in Ap- 
pendix [D|, to get, 

g[w] = const ^ wdvw (61) 
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We now make another change of variables from W a (x) to p(x), Z a (x), defined by equations 
(|?P and (H), to get, 

I = const J V[p]V[Z] \p,Z)(p,Z\ (62) 

where, 

V[p] = l[dp(x) V[Z] = Y[sin 2 9(x)sin0(x)d9(x)d4)(x)d4>(x) (63) 

X X 

Z has been parameterised as 

Z = (cos-e { 2 ),szn-e H 2 J) 

_ _ 

B. Overlaps 

The overlap of two gauge invariant coherent states \Wi) and IW2) is computed in the 
hydrodynamic limit in Appendix [E|. The final answer is, 



(Wi\W 2 ) = e 



-F[pi,Zi,p2,Z 2 l 



F = ~ l - I (pi(x) + p 2 (x))$(n 1 ,n 2 ) 



-7 / (PiW +P2(a:))hi( 1 



4Jx y ' y ' ' y " v 2 

-£(2n + 1) / ( Pl (x) - p 2 ( X ))( X \-L\y)( pi (y) - p 2 {y)) (64) 

where $(77,1,^2) is the solid angle subtended by the geodesic triangle with 774,712 and some 
third point on the unit sphere as vertices. 

Note that the overlap smoothly goes to 1 as (pi, Z\) — > (p 2 , Z 2 ). We will now evaluate 
the overlap for neighbouring states. We put p 2 (x) = p(x), pi(x) = p(x) + ed t p(x), Z 2a (x) = 
Z a (x) and Z la (x) = Z a {x) + ed t Z a (x). Keeping terms only to the order 0(e) we get, 

(W + ed t W\W) = e — 1p(^?(-)+o(^ 2 ) ( 65 ) 
where L 3 ^ = ^{Z^d^Z — h.c) for p, — t, 1 and 2. 
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If we now impose the LLL condition, then the charge density fluctuations get tied up to 
the the spin density fluctuations, i.e., 

The theory can then be expessed in terms of spin fluctuations alone. The expresion for the 
overlap in equation (^) then gets written as, 

(W + ed t W\ W) = e~ i€p $* L ^ x)+te ^+D L ^xLl(x)a v Ll(x)+o(e 2 ) ^ 

The second term in the exponent in the RHS of equation (^) is the Hopf term. Thus the 
theory, when restricted to the lowest Landau level is a NLSM with a Hopf term in the action. 

V. THE CLASSICAL LIMIT 

We will finally show that for large skyrmions, the theory becomes classical. Consider 
the set of states corresponding to configurations chaqracterised by a size parameter, A. We 
parameterie them as, 

Z\x) = Zfy (68) 

Substituting Pi(x), z^ a {x) and p^ix), Z2 a (x) in equation ( |64D and changing the variable 
x — > Xx we get, 

(W^W,) = eU x ^^W^ (69) 
\ ln( 1+n 2 1 '" 2 ) is zero h\ = fi2 and negative otherwise. Thus for W\ ^ W2, 

lim^lH^) (70) 

A^oo 

The coherent states thus become orthogonal when A — > 0. It can also be shown that the 
off-diagonal matrix elements of the observables in the coherent state basis, vanish in this 
limit. Hence the set of states corresponding to a system of skyrmions will behave classically 
in the limit of the skyrmion sizes tending to infinity. 



VI. CONCLUSION 



The motivation of this work was to examine the microscopic basis of the semiclassical 
NLSM for skyrmions at v — 1/ (2n + 1). In the systems of interest, the energy scales are such 
that it is important to impose the LLL constraint. The coherent state basis is the ideal one 
to address questions of quantum states in a semiclassical approach. Therefore, we developed 
a coherent state formalism of the composite boson theory. 

We showed that the coherent state basis of Hcb, when projected to the physical subspace 
Hphy, can be parameterized by a spinor field that we denoted by W a {x). In the hydrodynamic 
limit we have shown that these states, \W) themselves satisfy the coherent state properties 
of the resolution of unity and continuity of overlaps. The LLL condition is equivalent to the 
condition that W a {x) are analytic functions. 

The charge density is determined by the modulus of W i.e W^(x)W(x) and the spin 
density by the normalised CP\ spinor, Z a (x). In general these are independent quantities and 
therefore the charge density is independent of the spin density. However if W[x) is analytic, 
the modulus and phase of each of its components get tied up. We showed, that consequently, 
the excess charge density becomes proportional to the topological charge density which is 
determined by the spin density. Thus this proportionality will cease to hold in presence of 
Landau level mixing. We also showed that the condition for the total charge density to be 
proportional to the topological charge is weaker. It only requires W[x) to be analytic at 
infinity, i.e. that the ground state does not have Landau level mixing. 

Finally we showed that if we consider the set of states corresponding to classical con- 
figuations of characterised by a length scale A, then they become orthogonal in the limit of 
A — > oo. This implies that a system of skyrmions will behave classically in the limit of their 
sizes going to infinity. 
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APPENDIX A: THE PATH INTEGRAL REPRESENTATION 



In the following appendix we derive the path integral representation of the partition 
function by splitting the time interval t into N segments of length e and take the limit 
e — > 0, N — > oo such that eN = t. And at each intermediate step we insert the resolution of 
identity of H phy , 

Z = Tre^ Ht = Tr[e^ eH } N (Al) 

N N 

Z =U H(a n+1 <p n+1 \Pe=F m P\a n <p n ) (A2) 

n=0 ""'f"^" n=0 

where (a N+1 ,ip N+1 ) = (a , <p ) 

Since H commutes with P and P 2 = P, we get, after explicitly acting P on \a n (p n ) and 
making use of gauge invariance of \a<f)(a(p\, we get 

n . N 

Z=U „ IK (A3) 

where, 

t n = e ^l s ^ n ^(a n+1 - VP n+1 ^ n+1 e-^\e^ H \a n - VQ n - p n e-fe<°-+*>) 

(A4) 

To the order e: a n +i =«„ + ed„, Q n +i = + e^n and <^ n +i = + &Pn- And if we choose 
At+i — /^n + ^n+i then to the order 0(e) t n is: 

(a n - V/3 n - Vfi n + e(d n - VO n ), {y? n + e(0„ - £^ n fi n )} e -^ (/3 " +n " ) |e^ eH 
K - V/5 n - Vn n} ^e-t^+M) (A5) 

Using the fact that e/3 n = (3 n — fi n -\ + 0(e 2 ) and 

(a + 5a,(p + 5ip\e^ eH \a,(p) = (a + 5a,ip + 5(p\a,(p)[l - i^{oc(p\H\occp)] + 0(e 2 ) (A6) 
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where 5a ~ 0(e) and 5ip ~ 0(e) the above expression for t n , after defining a 0n = & n / c , 
and making use of gauge invariance of H, we get 



t n = exp[— J {--e^ vX 0i*{x)d v a*(x) + ea 0n (x)<^ n (x)v?„(x) - 

i%~[ip n (x)y n {x) - (f n (x)if n {x)] - H((f n , a n )}} (A7) 

If we now take the limit e — > the partition function becomes (after calling a by a) 

Z = J D[ao(x,t)]P[a,(x,t)]P[^(x,t)]e^ dM2a;£(a; ' t) (A8) 

where C(x,t) is the standard Chern-Simons lagrangian. 

APPENDIX B: THE WAVEFUNCTIONS 

In this appendix we give the details of the calculation of equation ([!(]). Using equations 
(H), (0) and (H we get, 

1 r r ^ 

lMK*<}) = — / C *i£/. a W xvn W(0|n[^)^(xi)I? t (^)]|a- Vn >V *-« n ) (Bl) 

Using the fact that, 

^(xOtffo) = e^+^-^i^O^^) (B2) 
and (|0|iC(a;) = 0, we can pull all the K's to the left and rewrite equation ( |BlD as, 

e-^/jSf=i a lW) x ( a >)- TO W) e 4^/j«>)-v^)-Er =1 «-W) 2 (B3) 

This is equation (|36l). We now write, 

<(x) = <' l J f,„U) (B4) 

where, 

- KX7 2 f m (x) = e5(x - x m ) (B5) 
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for m = 1 — N. 

The zero momentum mode of the Q integral will make the wavefunction vanish unless 
the total number of flux quanta equals the total numberof particles, i.e. 



k / V x a(x) = cN (B6) 

J X 

The VL integral is gaussian for the other modes and can be done exactly to give, 

1 K f N 

const x e xp[—- — / (oc(x) — ^«j(x)) 2 ] (B7) 
4 Tic Jx 

We then write the wavefunction as, 

M{Xi,<Ti}) = Const X e -|ll^)l 2 rj[ e ^-^)]2n+l 

i>j 

N N 

^[^(x^e-^^^le-^J-^-^x*^ 3 (B8) 



We also have, 



/ a(x).a m {x) = / Vnr(2c)-V/ m (ar) 

Jx Jx 

= - I n T (x)v 2 f m {x) 

Jx 



e -Q T (x m ) (B9) 



K 



Using the fact that the solution of equation ( p5|) is 



k 2ir 

and proceeding as in equation (|B9|), we have 



e 1 

f m (x) = — ln\x-x m \ (BIO) 



/ CX m (x).Oi n (x) fm(.Xr, 
Jx K 



e 2 1 



k 2vr 

So we have the result, 

N , e 2 1 Y 



o _ In \x m - x n \ (BH) 



/ (a(x) - £ ^{x)f = I vn T (x).vn T (x) - - 2 — E ln 

^ t=l ^ K Z7T m^n 



X-m X n 



-2- E ^(^m) + const (B12) 



e " 



K m=l 
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The (infinite) const comes from the m = n terms 
Finally the wavefunction is written as, 

^ N {{ Xi ,ai}) = const x V^)\\-hirJ x ™ T ( x yvn T ( x ) 

N 

x H( Zi - z,) 2n+1 Y[[^(x t )e^T(^n L{x ) }] (B13) 

i>j i=l 



APPENDIX C: THE OBSERVABLES 

We first evaluate the norm, J\f[W, Qt], in the hydrodynamic approximation, i.e., in the 
small amplitude and long wavelength limit. It is given by, 

N 

Af[w,n T ] = n / E\M{xi,°i})\ 2 

i=l Jx i Gi 

= const x e -/ I *.(*(*)e^ {nrW - ,lT(,), e - 1 fe/ I V!iTW'Vfl r W 

xflf e E l ln ^(^) H/ -(^)+ 2 >feE 1 ^(^)+(2r i +l)^ ij ln|x 1 -x j | 
i=l ' Xi 



As in the case of the Laughlin wave function, the norm has the form of a classical partition 
function of a 2D plasma. Here, there is also an "exernal potential" which is a function of 
W\x)W(x) and Qt(x). In the hydrodynamic limit, we write this partition function as a 
functional integral over the density field and evaluate it using the saddle point approxima- 
tion. So we change the variables from {:rj} — > p, where, 

~p{x)=j^5{x-x l ) (C2) 
i=i 

then for any function T of {x^}, we have 

N 



where the jacobian of the transformation is the entropy factor, 

J[p\ = e j>M-p(*) in pMI ( C4 ) 

Hence the norm can be written as, 
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' V[p]e- s[ P' w] (C5) 



where, 



S[p; W\ = [ [-p(x) + p(x) lnp(x) - p(x) \n{W a (x)W a (x)} - 2^-p(x)ft T (x) 



-2n(2n + l)p(x)^p(x)} (C6) 



We evaluate this functional integral in the saddle point limit. Dropping plnp term in 
comparision with p^p and substituting the solution of the saddle point equation, 

~p(x)=p- 47r(2 ^ + 1) V 2 ln{^(x)^(x)} (C7) 

we get equation (f46|). 

We evaluate the values of density and spin density by proceeding along similar calcula- 
tional steps employed in evaluating the norm, the density is given by, 

p(x) = (W\p(x)\W) 

1 N ( f \ - N 

= jy^F II ( J dx iH J 1pN({Xi,(Ti})^25(x - Xi^NdXi,^}) 

= i / V[p]p(x)e-^ 

= An(2~n + l) V2HWAx)Wa{x)} + f) (C8) 

where Z = J V[p] e - S ^ 
The spin density is, 

s a (x) = (W\s a (x)\W) 

1 N ( r \ - N 1 



J\f(W, 
\WAx)T aa ,WAx) 



2 W a n{x)W a n{x 

The current density is, 
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(p(x)) (C9) 



Jfa) = {W\J 3 {x)\W) 

1 N ( f \ - N 1 e 

+ /i.C 

= - -a*(z))] (CIO) 

c 

where i is the particle index and j is the coordinate index, a is defined through the relation 
kV x <3(x) = ep(x) and L| = ^(Z^djZ — h.c) for j = 1 and 2. 
The kinetic energy density is, 

T(x) = (W\f(x)\W) 

n(/^E)ElA^({x l ,a l })| 2 



1 = 1 \ CTj / j=l 



^ c {P(X)) (CU) 



where A = d*. + ~Zj and 2; = . 



APPENDIX D: THE T INTEGRAL 

In this appendix we evaluate Q[W] by doing the fix integral in equation ( T60|) . The 
saddle-point approximation of the integral gives 

Q[W] = const x e~L t^^f[W, fi T ] (Dl) 

where Qt is the solution to the saddle-point equation which, in long wavelength limit 
(V 2 \nWW < In WW) is, 

he 

Q T (x) = Q T (x) + — \n{W a (x)W a (x)} (D2) 
When this value for Qt is substituted in the above equation we get equation (|61~1). 

APPENDIX E: THE OVERLAPS 

The overlap of two gauge invariant coherent states \ Wi) and | W2), obtained by proceeding 
with steps similar to those involved in evaluating the norm, is 
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JAf(w 1 ,n T1 )Af(w 2 ,n T2 ) t\ \ J oj 



where 

f ab (x) = \n{W aa (x)W ba (x)} + 2 ]r- c Mx) (E2) 
If we express W in terms of p and Z we get the overlap to be 

(Wi\W 2 ) = e "M^)^ ln{Zlt(x)Z2(:!;)}v2ln{Zlt(a;)Z2(x)} 
x &\ Im {pl {x)+p2 {x)) ln{z i {x)Z2 {x)} 

The first term in the exponent of the RHS can be dropped with respect to the second one in 
the long wavelength limit. Making use of the relation, ZxaZ 2a = e^" 1 '" 2 ) ( 1+n ^ n ^ where 
$(ni, n 2 ) is the area of the spherical triangle with vertices at fix, n 2 and a third point on the 
unit sphere, in the above equation we get the equation (|64j). 
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